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Abstract 



We recently found that marine bacteria Vibrio alginolyticus execute a cyclic 3-step (run- 
reverse-flick) motility pattern that is distinctively different from the 2-step (run-tumble) 
pattern of Escherichia coli. How this novel swimming pattern is regulated by cells of V. 
alginolyticus is not currently known, but its significance for bacterial chemotaxis is self- 
evident and will be delineated herein. Using an approach introduced by de Gennes, we 
calculated the migration speed of a cell executing the 3-step pattern in a linear chemical 
gradient, and found that a biphasic chemotactic response arises naturally The implication 
of such a response for the cells to adapt to ocean environments and its possible connection 
to E. coWs response are also discussed. 

Keywords: V. alginolyticus, chemotactic behavior, biphasic response, motility pattern, drift 
velocity, chemotactic coefficient 
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Existing observations made in Escherichia coli (E. coli) have shown that sensing and 
motility impose different requirements on bacterial chemotactic response (1, 2). The debate 
on this interesting issue was initiated by the observation of Block, Segall, and Berg (3) who 
discovered that the experimentally measured chemotactic response function R(t) integrated 
over time t is zero. In physical terms R(t) can be thought as the Green's function of the 
chemotactic network when subjected to an impulsive or a 5-in-time perturbation. The im- 
portance of this null integrated effect goes without saying, and was immediately recognized 
by the investigators as the bacterium's means of sensing. In their words (3), "the bacterium 
compares the information received in the past one second with that received over the pre- 
vious three seconds." In effect, the double-lobe response function, which is displayed in 
Fig. d{c), allows the bacterium to react to fast temporal variations of a chemical signal 
c(t) but not to its dc component, enabling the cell to adapt to a wide range of chemical 
concentrations. Using a macroscopic diffusion argument, it was suggested by Schnitzer et 
al. (4) that a finite memory time is required for a bacterium to migrate in a linear chemical 
gradient; without the memory effect (or R(t) ~ S(t)), it was concluded that the chemotactic 
coefficient k = V/Vc or the drift velocity V would be zero, where V is in the direction of 
the chemical gradient Vc. However, de Gennes pointed out that the macroscopic diffusion 
approach ignored important correlations between bacterial swimming and the underlying 
chemical gradient (5). By taking into account such correlations, de Gennes showed that the 
optimal (or a fast) response for migration in a linear gradient is an exponential function 
with a decay rate determined by the cell's memory time r. He further pointed out that the 
double-lobe response function observed in E. coli could only reduce the migration speed in 
the gradient. 

Contributing to this stimulating debate is the finding of Clark and Grant (1), who argued 
that while a cell needs a fast drift speed in a concentration gradient, it is equally important 
for the cell to localize once the top of the gradient is reached. They showed that the single- 
lobe function proposed by de Gennes is inadequate for cell localization. By imposing the 
co-requirements of being able to localize as well as to migrate, they demonstrated that 
the optimal response function is biphasic, which is in remarkably good agreement with the 
one measured in the experiment (3). This observation led Clark and Grant to conclude 
that the biphasic response in E. coli perhaps reflects a compromised need of the cells in 
different environments. A recent study also suggested that the laboratory observed bacterial 
response corresponds to the maximin strategy that ensures the highest minimum uptake of 
nutrient for any profile of concentration (2). 

Recently, we found that the swimming pattern of the marine bacterium V. alginolyticus 
is a cyclic 3-step process (8), where a cell swims forward for a time interval Af and it then 
backtracks by reversing the motor direction for a time A&. However, upon resuming forward 
swimming, the bacterial flagellum flicks causing the cell body to veer in a new direction. 
This type of motility pattern is very different from that of E. coli, which exhibit a run- 
tumble pattern. By way of introduction a typical trajectory of V. alginolyticus and of E. 
coli are presented respectively in Fig. Ufa) and (b). For the V. alginolyticus' trajectory, Fig. 
QJa), the forward and the backward segments are designated by green and red, respectively, 
for clarity. We termed this novel swimming pattern a run-reverse-flick process. The last 
(flicking) step is functionally equivalent to a tumble in E. coli, allowing the bacterium 
to randomly select a direction, and a new cycle ensues. Despite the fact that run and 
reverse intervals, Ay and A&, as well as the flicking angle A9 are stochastic, the 3-step 
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Figure 1: Bacterial swimming trajectories of V. alginolyticus (a) and E. coli (b). The cells 
have been selected amongst many because they are more or less swimming in the focal 
plane, 50-100 //m above the glass coverslip. The starting points for both trajectories are 
indicated by the two large solid dots. The time lapse between adjacent dots are 0.067 s 
and 0.13 s for (a) and (b). The red and green segments in (a) designate the backward 
and forward swimming intervals, and transitions from backward to forward cause flicking, 
randomizing the swimming direction. Unlike a transition from forward to backward, which 
has a directional change A9 ~ ir (or backtracking), a backward to forward transition is 
random with A9 uniformly distributed between and 180°. (c) A hypothetical response 
function R(t) of E. coli based on Tu et al.'s model (6) is plotted as the dark curve, and 
the cell is stimulated at t = with a brief pulse of attractant (red curve). Here, R(t) = 



R 



exp 



I'm V Tf-fi 

methylation time (r„ 



■ exp 



where we set Rq = 1, and used the typical E. coli 



3 s) and phosphorylation time (r z = 0.5 s) (6, 7). 
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cycle is deterministic and has been observed in different V. alginolyticus strains and in a 
swimming buffer with and without a chemical gradient (8). In a steady state without a 
chemical gradient, we found that the probability density functions P(At) and P(Af) are 
statistically independent and have long exponential tails (or a Poissonian-like behavior) 
with the mean intervals r& — 77 ~ 0.3 s. However, when a point source of chemoattractant 
is present, the cells can quickly migrate along the gradient and form a tight pack around 
the source. 

The biochemical network that regulates the activity of E. coli motor is reasonably well 
understood (9). While this is not the case for V. alginolyticus, it cannot deter our progress 
because we know that even for very diverse microorganisms, such as E. coli and Bacillus 
subtilis that are roughly one billion years apart according to the recently constructed 
phylogenetic tree (10), the fundamental mechanism of regulation is still similar, i.e. a 
ligand binding to a receptor triggers a cascade of chemical reactions. The end product of 
the reaction is a chemically modified protein, called the response regulator (CheY-P), that 
binds to the motor, causing it either to rotate CCW (B. subtilis) or CW (E. coli). The 
basic aim of different microorganisms is also the same, namely guided by chemical signals, 
the cell is directed towards the source of chemoattractant and away from chemorepellent. 
According to the phylogenetic tree (10), V. alginolyticus appears to be much closer to E. 
coli than B. subtilis, suggesting that there is much in common between these different 
bacterial species. Indeed in V. alginolyticus, one can identify chemotaxis genes that are 
largely homologous to E. coli with the exception of cheV that is absent in E. coli but is 
present in B. subtilis. A recent study moreover showed that the phosphorylated CheY in 
V. alginolyticus causes the polar flagellar motor to reverse the direction from CCW to CW, 
similar to E. coli (11). 

It is clear that the 3-step swimming pattern is significantly different from the well- 
studied 2-step swimming pattern of run and tumble, and it has strong implications for 
bacterial chemotaxis, which can be characterized by an effective diffusion coefficient D and 
a drift velocity V in the presence or absence of a chemical gradient. The calculation below 
illustrates that cells executing the 3-step swimming pattern can exhibit rich chemotactic 
behaviors, and the variations can be acted on by natural selection so that a particular 
response emerges. Below we will illustrate these new aspects of bacterial chemotaxis based 
on our findings of the 3-step process. 

Similar to cells of E. coli, the flagellar motor of V. alginolyticus has two lifetimes for 
the state of rotations: one (r/) for the CCW interval and one (t&) for the CW interval, 
where the subscripts / and b stand for forward and backward swimming, respectively. To 
modulate their chemotactic behaviors, these lifetimes are affected by the local concentration 
of chemoeffectors and cells' adaptation mechanism. Unlike E. coli, however, CW rotation 
in V. alginolyticus causes the cell to backtrack. Both swimming intervals are expected to 
depend on the ligand concentration c(t), which we assume to be chemoattractant. For small 
c(t), we assume that a linear response is applicable and hence, 
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where Tf and tj, are the steady-state values, and Rf(t) and Rb(t) are the memory (or the 
response) functions, which are not necessarily the same for the two swimming intervals. 
In the above, an exposure to the ligand causes the forward lifetime to increase, and is 
consistent with our observations in V. alginolyticus (8). Linearity of Eqs. Q] and [2] suggests 
that it is possible to examine one delay time at a time and sum up all possible delays 
at the end. Following de Gennes, we write R s (t) = a s 5(t — 8), where the strength of the 
response a s (s = /, b) has the dimension of volume. Next, we consider a cell moving in a 
chemical gradient as depicted in Fig. [2j Our aim is to calculate the displacement Xj along 
the gradient in one cycle, Af + Ab, which leads to a mean drift velocity V = x~i/(' 7 7+ T fr) after 
averaging over Af and A&. Because a cell randomizes its swimming direction at the end 
of the backward interval by a flick, the motions in two consecutive cycles are uncorrelated. 
This allows us to place the origin of time (t = 0) at the beginning of the forward run. 
Assuming that the forward run time is Poisson distributed, the surviving probability of a 
cell swimming forward up to Af is given by, 



P f (A f ) = exp 
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1 



T f {t>) 



exp( 1 - i 



i + / dt ' c {t') 



(3) 



and the probability that it stops immediately after Af is —dPf(Af)/dAf. Likewise, the 
surviving probability of a cell swimming backwards from A f to A f + A& is given by, 



P b {A b ,A f ) =exp 



A f +A b 



dt' 
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- exp( ) 
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and the stopping probability at the end of the backward run is — dPb(Ab, Af)/dAb- 
follows that the net mean displacement in one cycle is given by, 



(4) 



It 
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where Xfi and x b i represent respectively the mean displacement during the forward (Af) 
and the backward (A^) swimming interval, and (...) designates the angular average for 
Vfi and Vbi- For the linear gradient depicted in Fig. [21 the concentration experienced 
by the cell can be represented as c(t) = cq + Vc ■ Vfi ■ t for < t < Af and c(t) = 
Co + Vc ■ Vfi - Af + Vc • • (t — Af) for Af < t < Af + Af,. Since cq is determined by the 
velocity in the previous cycle, it does not contribute to the above integrations after angular 
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averaging. Although the calculation of Eq. [5] is tedious, which is given in Appendix A, the 
final result is straightforward: 
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The first term in the curly bracket of Eq. [6] is the displacement during the for- 
ward interval, and the second term is the displacement during the backward interval. 
It is noteworthy that during the second interval, there is a cross term proportional to 

T f-T b ( v fi v bi) ex P( — Tf) ~ 7J ex P( — Tfc))> wri ich results from the delay, i.e., even though 
the cell is moving backwards, in the early episode of that interval, the cell still remem- 
bers the concentration sensed during the previous forward swimming. This gives rise to 
anti-correlation, since {vfiVbi) < 0, that contributes to a negative displacement. This 
important correlated motion adds richness to bacterial chemotaxis and is what makes 
V. alginolyticus behave differently from E. coli. We noted that Eq. [6] yields the result 
xi = a/r|(ujj) + afer;, (t&(u^) + Tf(vfiVbi)) Vc in the limit of no memory, 9 — > 0. It is 
interesting that even when there is no memory, the cross term survives because there is no 
direction randomization after a forward run. Moreover, the total displacement during the 
backward interval, can contribute positively or negatively to the displacement, depending 
on the mean lifetimes Tf and T b , and the swimming speeds Vf and Vf,. We found the swim- 
ming pattern of V. alginolyticus is approximately symmetric with Tf ~ Tb and \vf\ — \vb\ 
(8), and hence X{ ~ afT 2 (v'j i )'Vc. 

The corresponding quantity for E. coli is x~i = afTf(v'j i )(/)oVc when MO, where the 
subscript / stands for the forward run (or CCW rotation) and <pQ = tqcw I ( T CCW + T cw) 
is the CCW bias. Since near a steady state 4>o ~ 0.5 or ~ 0.8 according to Refs. (3) and 
(12), respectively, E. coli cells produce a smaller displacement than V. alginolyticus within 
one swimming cycle if everything else is equal. Using Vfi = —Vbi = Vi and summing up all 
possible delays, we found from Eq. [6] that the mean displacement is given by 
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The average drift speed in the gradient is V{= 
chemotaxis coefficient k to be calculated. In E. 
coefficient D 



kS/c) = Xij{jf + Tb), which allows the 
coli, k is proportional to the diffusion 
^(j)o(v 2 )Tf and one finds k = D Rf(9) ex.p(—^j)d9. For an organism 



exhibiting the 3-step swimming pattern, the diffusivity is given by: 
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Figure 2: Migration of V. alginolyticus in a linear chemical gradient, (a) In the spatial 
domain, the chemical gradient is specified by the dashed blue line. The green and the red 
arrows correspond to forward and backward swimming segments along the gradient. The 
color usage in the arrows is consistent with that in Fig. U Aj and are respectively the 
forward and backward swimming time intervals, and Vfi and vu are respectively the forward 
and backward velocity components along the chemical gradient. Note that backtracking 
means v\, = —vt. (b) The bacterial chemotactic network processes the chemical information 
in the temporal domain, and the concentration detected by the cell is depicted in (b), where 
9 is the memory time of the bacterium. I and II are chemosensing in the current cycle, and 
I' is due to the previous cycle. 
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and the chemotaxis coefficient can be written as 
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This calculation leads to two possible scenarios (or fundamental hypotheses) for bacterial 
chemotaxis: (i) independent and (ii) shared chemosensing. In the first case, the response 
functions in the forward and backward intervals are uncorrelated, i.e., Rf(9) and Rb{9) have 
different functional forms, so that the sensing system breaks the time reversal symmetry. 
In order to achieve such a control, the flagellar motor cannot only passively receive signals 
from the chemotaxis network but instead the status of the motor must be made known 
to the chemotaxis regulatory network. This may be attained either by the flagellar motor 
being a part of the regulatory network or by a feedback signal via a protein that can reset 
the chemotactic response. In short, there will be a back flow of information from the motor 
to the chemotaxis network in addition to the normal chemotaxis regulation. To optimize 
the drifting velocity, we applied a variational principle to Eq. [9J which is delineated in 
Appendix B. We used the constraints that Rf(9) and Rb(9) have constant variances ct1/t s 
{s = /, b) (1), yielding- 
ly^) oc ^-exp( ), (10) 



R b (0) oc ^ 
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exp( ) — — exp( ) exp( 1 
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It is evident from the optimization procedure that in order to attain the maximum 
possible drifting speed, the forward response function Rf(9) should be monophasic but 
the backward response function Rb(6) can be either monophasic or biphasic, depending on 
the ratio of the two lifetimes, f3 = Tb/T f . Fig. [3] displays (r/ /af)Rf{9) and (T b /ab)Rb(9) 
for different values of (3 = 0.8, 1.2, 1.5, 1.8, and 2.4. The figure shows that the biphasic 
character of Rb(0) becomes more pronounced as /3 decreases towards unity but dissappears 
altogether for j3 < 1, where the response is negative for all 9. An analysis shows that the 
biphasic response occurs in a narrow range of f3 (1 < ft < 2), and outside this range the 
response is always monophasic. This behavior is understandable since when Tb is shorter 
than Tf, the backward interval is strongly influenced by the signal sensed in the previous 
forward interval due to the memory effect. To deal with this inconsistency between sensing 
and motility, the optimal strategy is a negative monophasic response as depicted by the 
purple curve (/3 = 0.8) in Fig. [3l On the other hand, when is longer than Tf, the 
cell would have consistent sensing and motility so that a monophasic positive response is 
more favorable, which is shown by the red curve (/3 = 2.4) in Fig. [3l In the limiting case 
Tb S> Tf or Tb Tf, Eqs. [9] QUI and [TT] make it clear that the chemotactic coefficient k is 
dominated respectively by the backward or the forward swimming interval. The situation 
is formally equivalent to E. coli chemotaxis, where the monophasic response is optimal for 
a fast migration in a linear chemical gradient as was concluded by de Gennes (5). 
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Figure 3: Chemotactic Strategy I. The bacterium uses separate response functions, Rf(9) 
and Rb(9), for chemosensing. The figure shows the dimensionless forms of the response 
functions. Here the black curve is for Rf(0) and the purple, blue, pink, green, and red 
curves are for Rb(0) with /3(= t^/tj) = 0.8, 1.2, 1.5, 1.8, and 2.4, respectively. The inset 
is the phase diagram for chemotactic strategy II. It displays the phase boundaries between 
monophasic (I, III, IV) and biphasic (II) response regimes when the chemotaxis response 
obeys the relation R(6) = R f (8) = Rb{0)h- 
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In the case of shared chemosensing, the bacterium uses a single response function R(9), 
albeit the amplitudes of the responses may be different in the two directions, Rf{6) = 
Rb{9)/l = R(ff). A simple reason for 7 ^ 1 could be due to different swimming speeds Vf 
and Vf,, but other possibilities may also exist. For this type of sensing, there is no breaking 
of time reversal symmetry since the chemotaxis network processes information received 
during the forward and the backward interval equally, and there is no need for a back flow 
of information. Using Eq. [9j we found: 
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Applying the variational principle again (see Appendix B), we found that the drift velocity 
is optimized by the following response function 
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As displayed in Fig. R(9) can be monophasic or biphasic depending on 7 as well 



as the time ratio ft = 17, /ry. The biphasic regime is bounded by 



1-/3 
P 



< 7 < 
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for 



< f3 < 1 and 7 > for 1 < (3 < 2, which is displayed in the inset of Fig. [3l The inset 
shows that the parameter space (7, /3) consists of four different regimes with I, III, and IV 
being monophasic and II biphasic. Our theory hence predicts that if a bacterium uses a 
single response function, for very short (/3 <C 1) or very long {(3 3> 1) backward swimming 
intervals, the biphasic response is not a good chemotactic strategy for migration in a linear 
chemical gradient. The biphasic response emerges only when Tf and T& being close (or 
P ~ 1), which is the case in V. alginolyticus (8). It is conspicuous that in the limits /3 — >■ 1 
and 7 — > 1, R(9) calculated using Eq. [13] is identical to the solution of a critically damped 
harmonic oscillator, which has the interesting property of Jq°°(1 — 9 /r) exp(— # /t)(19 = 0, 
i.e., the response is "precisely" adaptive. 

The above two hypotheses are testable by laboratory experiments where the bacteria 
are subject to a defined chemical stimulation, and one measures the switching rate s(t) = 
2/( r /(0 + T b(t)) and the forward swimming bias (j){t) = Tf(t)/(Tf(t) + Tb(t)) as a function 
of time. For a weak stimulation, the above calculation allows us to find, 
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Figure 4: Chemotactic Strategy II. The bacterium shared the same response function 
R(9) = Rf(9) = Rb(0)/^f for the forward and backward swimming intervals. In (a), 
P = 0.75, and 7 = 0.1 (red), 1.0 (pink), 2.0 (blue), and 3.0 (purple). In (b), 7 = 2, 
and j3 = 0.1 (red), 0.3 (pink), 0.5 (blue), 1.0 (purple), and 2.0 (green). As can be seen, 
for fixed (3 ~ 0.75, the response becomes strongly biphasic as 7 increases. On the other 
hand, for fixed 7 = 2, the response is monophasic for small /3, and becomes biphasic for 
intermediate values of /3, and returns to monophasic for /3 > 2. 
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where so = 2/(rf + T b ) and (j>o = Tf/{rf + T b ) are the steady-state switching rate and the 
forward bias, respectively. The expressions are significantly simplified if the perturbation 
is (5-in-time, c(t) = c'5(t), and they are given by, 

s(t) = s [1 - c' (<p Rf(t) + (1 - <Po)R b (t))] , (16) 

4>{t) = ^[l + c'(l - <fo) (R f (t) - R b (t))) . (17) 

The calculation shows that if the second scenario is true and Rf(t) ~ R b (t), the forward 
bias will be weakly dependent on the time t, and the switching rate is simply given by 
s(t) ~ So[l — c'Rf(t)]. However, if the first scenario is true, the measured s(t) and 4>(t) can 
be used to find the response function Rf(t) and Rb(t) using Eqs. [T6l and [T71 In this case, 
the following simple relations result, 



R b (t) = ~~7 



<m s(t) 

l - <t>(t) s(t) 



(18) 
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An alternative and perhaps more direct way to find Rf(t) and R b (t) is to perform 
conditional stimulation for individual cells. The bacterium can be either tethered to a 
surface, such as in Block et al.'s experiment (3), or freely swimming, as in Khan et al.'s 
experiment (13). For tethered cells, one can apply a pulse of chemoattractant at the moment 
the motor switches from CW (CCW) to CCW (CW), and record the subsequent swimming 
interval A\f (A^), where the subscript 1 emphasizes the interval before the first switch. By 
counting the switching events up to time t, one can construct a cumulative PDF (normalized 
by the total number of cells) *& s (t), and the time-dependent switching rate can be obtained 
according to T~^(t) = — ^ln(l — ^f s (t)), where s = f or b. For freely swimming cells, 
one can use photo-active serine, which is an attractant to V. alginolyticus, to stimulate 
cells. If the first scenario is true, one should find that tJ (t) and T h ~ l {t) have different time 
dependence or equivalently Rj(t) and Rb(t) have different functional forms. However, if 
the second scenario is true, there should be not much difference between Tj l (t) and T fe ~ 1 (t) 
or R f (t) oc R b (t). 

To conclude, the 3-step motility pattern of V. alginolyticus discovered in our recent 
experiment (8) allows significant variations in bacterial chemotactic behavior. These varia- 
tions can be acted on by natural selection and give rise to distinct phenotypes observed in 
the wild. Compared to the 2-step swimming pattern of E. coli, cells of V. alginolyticus can 
engage in chemosensing and migration in both the forward and the backward swimming 
intervals, and hence their "duty cycle" is ~ 100% as compared to ~ 50 — 80% in E. coli (3, 
12). An important aspect in 3-step chemotaxis is backtracking that gives those bacteria 
heading down a gradient an opportunity to re-exploit what they find a moment earlier. In 
our opinion, the full duty cycle, backtracking, and flicking are defining characteristics of 
V. alginolyticus. These significant niches are likely selected for by the ocean environment 
where a quick response to transitory signals is important. We showed that for a swimmer 
executing the cyclic 3-step motility pattern, a biphasic response arises naturally without 
the need to invoke cell localization as suggested for E. coli (1). Moreover, we showed that 
the biphasic response is most effective when the forward rt and the backward r b swim- 
ming intervals are comparable. This makes biological sense since a brief forward or a brief 
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backward interval contributes little to motility, and consequently a monophasic response is 
sufficient for migration. This also raises the interesting question why the non-motile CW 
interval in E. coli is so long, taking up at least 20% of the duty cycle. If tumbling is just 
to change the direction, would not it be better if CW interval is shorter? An interesting 
possibility is that the ancestral cell that gave "birth" to E. coli and V. alginolyticus was a 
3-step swimmer. However, when E. coli became specilized in a different environment, which 
favored multiple flagella for motility, they gave up backtracking and flicking, resulting in a 
tumbly movement. In this view, then, it is not surprising that E. coli's tumbling interval 
is long and its chemotactic response is biphasic. 

Based on motility alone, we propose two different mechanisms, independent and shared 
chemosensing, by which cells of V. alginolyticus can optimize their migration speed in a 
linear gradient. Interestingly, the biphasic response appears in both types of chemotactic 
strategies. We know very little at present how the chemotaxis network of V. alginolyticus 
regulates the 3-step motility pattern. However, based on our calculation, it is likely that 
a biphasic response is also adopted by V. alginolyticus, and it is awaiting to be verified 
in future experiments. Finally, it would be interesting to generalize the above calculation 
to situations where the chemical landscape is constantly changing, such as chemical waves 
(14). 

This work is supported by the NSF under the grant no. dmr-bp0646573. 



Appendix A: Calculation of the Mean Displacement 

In the following we provide a more detailed derivation of the mean displacement X{ = 
Xfi + in a single 3-step cycle. The displacement is made in the two time intervals, 
Af and and is represented by Eq. G3 The concentration sensed by the bacterium is 
piecewise continuous according to Fig. [2] and is given by, 

t < 

< t < A f 
A f <t< A f + A b 

where the subscript i designates the component of the velocity along the gradient direction. 
The primed and unprimed velocities correspond to t < (regime I') and t > (regimes I 
and II), respectively. 

The first part of Eq. [5] is readily calculated by integration by parts, 

*/*=(/ dA f (-^r) vfAf ) = \j dA f p f( A f>fi) ( 2 °) 

^0 ' \o ' 

where Pf(Af) is given by Eq. [3j which contains an integration in time t over the range: 
—8 < t < Af — 8. Since Af varies from to oo, we have to distinguish two cases in the 
integration: (i) Af — < and (ii) Af — 9 > 0. One can deal with these two mutually 
exclusive cases by the use of Heaviside functions H(x), i.e., we write, 



c{t) 



c + Vcv' u t , 
c + Vcvfit, 

c + VcvfiAf + Vcv bi (t - A f ) , 
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oo 



I dA f (...) = J dA f [H(0 - A f ) + H(Af - 9)}(...). (21) 
o o 

The first Heaviside function confines the integral to t < 0, and since {v'^Vfi) = 0, there is no 
contribution from this term. The integration constrained by the second Heaviside function 
yields, 



x fi 



a f Vc(v 2 fi )T 2 f exp(-^-y (22) 

This equation is identical to that found by de Gennes when he calculated the drift velocity 
for E. coli cells (5). 

The second part of Eq. [5] is more complicated because one has to take into account 
more possibilities. Again, we used integration by parts to obtain, 

*. s ( ,, M , M -^)(-^)^\ 



o 

oo 



= (p f (0) j dA b P b (A b ,0)v bi J + (^j dA f j dA b P f (A f ) dPb{ ^ Af) v bl j . (23) 

Let the first term in the above equation be xffl and the second term be xj^ '. Since 
Pj(0) = 1 and P b (A b ,0) = Pf(A b ), it follows that the integration in the first term is 
identical to Eq. [20] with the replacement of the subscript / by b. This yields, 

4 b) =a 6 V C (4)r 6 2 exp(-^. (24) 

Now, lets examine the anti-correlation term , which corresponds to the situation when 
the bacterium swims down the gradient but it still keeps its "old good memory". Dropping 
the nonlinear terms in concentration c, we found, 

, oo oo Af+A b —9 

4 f) = y(f dA f J dA b exp (-^j exp (-^j JL J dtc(t) Vbl \ . (25) 

When integrating over Af, there are two possibilities for the lower limit of the ^-integration, 
i.e. either Af — 9 < or Af — 6 > 0. These will be delimited by the Heaviside functions as 
before. For each of these cases, while integrating over A b , there are additional possibilities 
for the upper limit of the i-integration. For the first case, when Af — 9 < 0, there are 
three possibilities: (i) Af - 9 < A f + A 6 - 6 < 0, (ii) 0<A f + Afe-6><A/, and (iii) 
Af < Af+A b — 9 < Af+A bl corresponding to the regimes I', I, and II in Fig. [2j respectively. 
However, since motion is uncorrelated after a flick or (vi-Vf^) = 0, the first possibility does 
not contribute to the displacement. In the second case, when Af — 9 > 0, there are two 
additional possibilities: (iv) A f -8 < A f +A b -8 < A f and (v) A f < A f +A b -9 < A f +A b , 
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corresponding to the regimes I and II in Fig. respectively. The corresponding time 
integrals for the above four possibilities (ii-v) are given by 



dt c(t) 



c {A f + A b -9) + lVcv fi (Af + A b 
c A f + \Vcv fi A) + ci (A 6 - 9) 
+\Vcv bi {A b -9f , 
c A b + Vcv fi A b (A f + \A b - 9) , 
c 9 + Vcv fi 9 (A f - I) +ci (A 6 - 9) 
+\Vcv bi (A b -9f , 



, (ii) 

(in) 
(iv) 

(v) 



(26) 



where c\ = cq + VcvfiAf. Using the above expressions, we take the derivative with respect 



to A j to obtain, 



A f +A 6 







dA t 



dt c(t) 



Ar-e 



c + Vcv fi (Af + A b -9), (ii) 

c + Vcvfi(Af + A b - 9) , (Hi) 

Vcv fi A b , (iv) 

Vcv fi A b . (v) 



(27) 



Again, using the Heaviside function to represent these four non-trivial possibilities, we have 
the following identity 



dA f J dA b (.. 



dA f I dA b 



{H(9 - A f ) [H(A f + A b - 9)H(9 - A b ) + H(A b - 9)H(9)} 
+H(A f - 9) [H(A b )H(9 - A b ) + H(A b - 9)H(9)]} (...). (28) 



Substituting this equation into Eq. [25l we found, 



- r {bf) 



n 



oo oo 

dA f / dA b exp ( ) exp ( 

J V T i ) V n 





x {H(9 -A f ) [H(A f + A b - 9)H(9 - A b ){ii) + H(A b - 9)H(9)(iii)] 

+ H(A } - 9) [H(A b )H(9 - A b )(iv) + H(A b - 9)H(9)(v)}} v bi ), (29) 



where (ii), (Hi), (iv), and (v) are the terms given in Eq. [271 The four integrations in 
the above equation are delimited by different combinations of Heaviside functions, yielding 
different lower and upper integration limits for each integral. Designating these integrals 
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I u \ l III \ 



, and ( x^P ) , we found, 



V 

^-Vc(v fi v bi ) J dA f J dA b exp(-^jexp(-^\(A f + A b 

B-A f 



a b Vc{vfiV b i) 



Tf ~ T b 



■ cxp 



1 1 



X Ub 



exp(-)-l 



(30) 



-AM) 



e oo 



Vc(v fl v bi ) J dA f J dA b exp f-^J ex P (~~) ( A / + A &-#) 

9 ^ b 

r f + n-(T f + T b + 9)exp(^-^- S j exp^-^J, (31) 



a b S/c(v fi v bi )Tf 



oo e 

^-X7c(v fi v bi ) J dA f J dA b exp ex P ( _ ~^~) Afc 



a b Vc{vfiV bi )Tf 



r b - (9 + T b ) exp 



n 



cxp 



(32) 



oo oo 



;<&/} 



< y-Vc(v fi v bi ) J dAf J dA b exp ^-^-^j exp A 6 
a b Vc (vfiv bi ) Tf(n + 9) exp ex P (~^) ' 



(33) 



In the above calculation, the terms involving cq do not contribute since (i^) = 0. The 
anti-correlation term due to all the above contributions is then given by, 



=<V> 
l bi 



-M) 



i4 v, ).. + (4 w) 

/ in \ 



a b Vc {vfiV bi ) 



T 2 T 2 
T f T b 

n-Tf 



1 

V 



+ 


H w> ) 











(34) 



Combining Eqs. [22j [JU and [341 we finally obtain the mean displacement in a given cycle 
for the 3-step swimmer, 



a/Vc {vfyTf exp ( -^-J + a 6 Vc (v^) r fc 2 exp — 



r 2 r 2 
T f T b 

T b -Tf 



1 / 0\ 1 
— exp exp l 

T f V nj n V r / 



(35) 
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Appendix B: Drifting Velocity Optimization 

For the first chemotactic strategy, Rf(t) and Rb(t) in Eq. [9] are independently optimized. 
The procedure requires to constrain a family of response functions R s (t), where s = f, b. 
We followed Clark and Grant's approach (1) and assumed that R s (t) is finite, continuous, 
and decays to zero for large t. The simplest way to impose the constraint is to assume a 
finite variance 



R z s (t)dt = a;/r s (36) 
n 

that is to be satisfied by all curves in the family. Optimizing k with the above constraint 
is equivalent to 



SR s (t) 

o 



oo 

J dt R s (t)K s (t)-\(R 2 s (t)-^] 



(37) 



where K s {t) is the kernel that weights the forward (s = /) and the backward (s = b) 
response functions, 

K f (t)=e X p(-^ r y (38) 

K b (t) = exp(--) - — T J— (I e xp(--) --exp(--)) . (39) 

Aside from normalization constants, the optimized response functions are given in Eqs. [TU1 
anddB 

The similar procedure can also be applied to the second chemotactic strategy, resulting 
in the optimized response function given by Eq. [T3l 
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